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Gradient index devices for the simultaneous control of the So and Ao Lamb modes are presented 
and numerically analysed. It is found a special relationship between the refractive indexes for 
these two modes which allows their simultaneous control by means of graded phononic crystals and 
thickness variations of the plate. A flat gradient index lens and a circular Luneburg lens are designed. 
Numerical simulations show that the performance of these devices is good for the two modes in a 
broadband frequency region, and that this approach can be used to design more advanced refractive 
devices for the total control of guided vibrational modes. 


Gradient index (GRIN) devices consist of locally inho¬ 
mogeneous materials in which the refractive index is a 
function of the spatial coordinates, therefore waves fol¬ 
low curved trajectories that can be properly designed to 
focus the energy at a specific point. These devices, which 
are well known in optics [I], are recently receiving increas¬ 
ing attention in the field of phononic and sonic crystals, 
given that these structures can be properly designed to 
work as GRIN devices. 

Effectively, phononic crystals are essentially periodic 
arrangements of solid inclusions in an elastic matrix (fluid 
background for sonic crystals) and they behave, in the 
low frequency limit, as homogeneous materials with ef¬ 
fective constitutive parameters which depend mainly on 
the filling fraction of the inclusions [2]. Then, by lo¬ 
cally changing this parameter according to some specific 
law, we can create an artificially inhomogeneous material 
which will work as a gradient index device. Following this 
approach, a wide variety of inhomogeneous devices have 
been proposed, like GRIN flat lenses EH2, acoustic beam 
modifiers [8], omnidirectional refractive devices [3 HU] or 
absorbers mm and acoustic or elastic cloaks p~3l IT4] , 
which are indeed a new type of advanced gradient index 
device. 

In the case of GRIN devices for solids, the main ap¬ 
plications arise in the propagation of guided waves, like 
Lamb or Love waves in plates and substrates [15]. How¬ 
ever, these systems present several propagation modes 
and polarizations, being their dispersion relations differ¬ 
ent in general. The main effort so far in this domain has 
been focused on the design of GRIN devices for the con¬ 
trol of the lower order Lamb mode in thin plates, i.e., for 
the control of flexural waves [TOUTS] . 

In this work, the performance of this type of devices 
is improved, by exploiting an interesting property of the 
two lower order Lamb modes, which is the fact that they 
are connected by means of the so called “plate velocity”, 
and it is shown that it is therefore possible to design of 
GRIN devices for the simultaneous control of these two 
modes in a broadband frequency region. 


The dispersion relation of Lamb waves in elastic plates 
can be found in reference [15j where it is also shown that 
in the low frequency limit, defined by means of wave¬ 
lengths A > h, being h the thickness of the plate, the 
dispersion relation for the So mode is linear and given by 

w = v p k s , (1) 

while for the Ao mode it is parabolic and given by 

,2) 

being v p the so called “plate velocity”. These expres¬ 
sions show the interesting property that the dispersion 
relation of the A 0 mode depends on the plate’s thickness 
h , and it is directly related with the dispersion relation 
of the So mode by means of the plate velocity v p . In this 
work this relationship between the dispersion relation of 
the two modes will be exploited, and homogenized two- 
dimensional phononic crystals will be employed to that 
purpose. 

The above expressions are valid for the vibrations of a 
homogeneous and isotropic plate, and it must be figured 
out if they can also be applied to phononic crystal plates 
in the low frequency limit. Phononic crystal plates con¬ 
sist of periodic arrangements of inclusions in a plate, and 
they can be considered like “slides” of a two-dimensional 
phononic crystal consisting of periodic arrangements of 
cylinders. This two-dimensional phononic crystal will be¬ 
have, in the low frequency limit, like a homogeneous but 
anisotropic material, given that the properties in the xy 
plane are different from those in the z direction. Then, 
a “slide” of this phononic crystal will not vibrate like 
a homogeneous and isotropic plate, but like a homoge¬ 
neous and anisotropic plate. It is not clear therefore that 
expressions ([l]) and Q be valid in this case. 

The analysis of Lamb waves in anisotropic plates is 
complex, but in the specific case of tetragonal materials 
(triangular or square arrangement of inclusions), which 
is the case considered here, it can be demonstrated (see 
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supplementary material) that expressions 0 and (J2| are 
valid but the plate velocity must be given by 


p p \ C11C33 j 


( 3 ) 


where Cij stands for the ij element of the stiffness tensor 
of the plate. 

The propagation of waves in a medium is determined 
by the dispersion relation, but the refraction of waves 
when passing from one medium to another one is man¬ 
aged by the ratio between the wavenumbers in the two 
media, which defines the index of refraction. Then, it is 
easy to check that if a wave changes from propagation 
medium 1 to medium 2, the refractive index for the So 
mode will be 


n s = 


Vpi 

Vp2 : 


while for the A 0 mode will be 


n A = n s 


hi 

V 


( 4 ) 

( 5 ) 


where it has been assumed that the two media are two dif¬ 
ferent plates, with different materials and thickness. The 
above two expressions show that, although the disper¬ 
sion relation for the two modes is different, with different 
frequency-dependence, their refractive indexes are closely 
related, with the remarkable result that, if the thickness 
of the two media is chosen such that hi = ft, 2 ^ 5 , it will 
be found that the refractive index of the So and the Ao 
mode is the same. This interesting property will be ex¬ 
ploited here for the design of a refractive device identical 
for the two modes. 

Two devices are proposed in this work based on the 
combination of graded phononic crystals and thickness 
variations. The first one consists in a gradient index lens 
and the second one consists in a Luneburg lens. These 
devices have been already studied for flexural waves, and 
here it is proposed their design to work simultaneously 
for the two fundamental Lamb modes. 

A GRIN lens consists in a rectangular device where the 
refractive index is a function of the distance to the lens’ 
axis, being such a function given by[3] 


n(y) =n 0 sech(ay), (6) 


so that if the gradient of the index is along the y direction, 
a wave travelling along the x axis and arriving to the 
device is focused at a distance / = n/2a. Similarly, a 
Luneburg lens consist of a circular device in which the 
refractive index is a function of the distance to the center 
of the lens and given b y ID 

n(r) = y/2 - (r/i? L ) 2 , (7) 


where Rl is the radius of the lens. This device is an om¬ 
nidirectional device in which any plane wave arriving at 


its surface is focused on the diametrically opposed border 
of the lens. 

Phononic crystal based GRIN devices consist in the 
realization of the above variations of the refractive index 
by means of a variation in the physical properties of the 
inclusions. In the low frequency limit phononic crystals 
behave like homogeneous materials with some effective 
parameters which are function of the inclusions’ proper¬ 
ties and sizes, then a rectangular or circular region of a 
phononic crystal can be properly designed to satisfy the 
above position dependent refractive indexes. 

The design of this type of devices requires therefore of 
a proper effective medium theory which gives the effec¬ 
tive refractive index of the phononic crystal as a func¬ 
tion of the physical properties of the inclusions. In the 
present case, this theory is the recently developed ho¬ 
mogenization theory for flexural waves (Ao mode) which 
can be used to obtain the refractive index of the Aq mode 
m- In the framework of this theory the propagation of 
waves is described in terms of the rigidity D 5 of the plate 
and its mass density and if in the plate there is a 
phononic crystal, the theory gives the effective rigidity 
D * and mass density ft* relative to that of the back¬ 
ground Db and p&. In the supplementary material it is 
shown that the plate’s velocity in the phononic crystal 
region is related with this theory by 

_2 = E b 

Vp p* p b (l - vl) 

so that from the homogenization theory for flexural waves 
it is possible to deduce the plate’s velocity and, conse¬ 
quently, the refractive index not only for the A 0 mode 
but also for the So mode. 

The design procedure for the GRIN lens is therefore 
as follows. First, an inclusion located at R a = (x a ,y a ) 
is selected, the expected refractive index at that position 
is computed by means of equation © and by means of 
equation 0 it is obtained the corresponding plate ve¬ 
locity. Equation 0 relates the plate velocity with the 
properties of the inclusion, since the objective is to create 
the GRIN lens by changing the filling fraction, from this 
equation we solve for the expected filling fraction and, 
subsequently, for the inclusion’s radius. The procedure 
is repeated to all the inclusions in the plate so that the 
design of the full GRIN lens is finally completed. Notice 
that it has been employed the effective medium theory 
for flexural waves to obtain the effective parameters, but 
equation 0 for the design of the lens, what means that 
the GRIN designed works only for the So mode. 

The key point in this work is that it is possible to 
choose now a position dependent thickness of the plate 
such that h 2 (y) = hb/n(y ), being h 5 the thickness in 
the background (i.e. outside the lens region), given that 
n(y) = ns{y) and following equation ([ 5 ]), it will be found 
that the refractive index of the two modes is identical 
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and, therefore, the GRIN lens will work properly for the 
two modes. 

The designed lens consists of a square arrangement 
with lattice constant a of lead inclusions in a silicon plate 
of thickness h = a, with p a = 11.34 kg/m 3 , E a = 16 GPa 
and v a = 0.44 the mass density, Young’s modulus and 
Poisson’s ratio for lead and p 5 = 2.33 kg/m 3 , E & = 150 
GPa and = 0.28 the corresponding parameters for sil¬ 
icon. In the low frequency limit this arrangement of in¬ 
clusions will behave like a tetragonal material, but it has 
been mentioned before and demonstrated in the supple¬ 
mentary material that the basic relationships used here 
work as well for this type of plates. The rectangular re¬ 
gion of the phononic crystal is selected as having N y = 15 
rows and N x = 30 columns, which give us a lens of height 
L y = 15a and length L x = 30a. 

Figure [l] shows COMSOL simulations of the interac¬ 
tion of the GRIN lens with the So and Aq polarizations. 
Upper panels show the proper GRIN design in which the 
thickness of the plate in the rectangular phononic crys¬ 
tal region is position-dependent, while mid panels show 
the GRIN lens without the mentioned thickness varia¬ 
tion. Left panels correspond to the the normalized z 
component of the displacement field for the A 0 , and the 
right ones to the So mode, in which it is plotted the nor¬ 
malized x component of the displacement field. Given 
that the performance of the device depends on the wave¬ 
length of the field and not on the frequency, simulations 
are shown for frequencies in which the two modes have 
the same wavenumber ka = 1, which corresponds to 
uja/v p = 1/vT2 for the Ao mode and uoa/v p = 1 for 
the So one. As can be seen, the device works properly 
for the two modes, being the focusing point identical in 
both situations. It can also be seen that the field dis¬ 
tribution is not the same for the two modes, given that, 
although the Ao mode can be scattered almost without 
interaction with other Lamb modes, the So mode tends to 
excite SHo waves, hence the displacement field is not the 
same. In other words, this theory provides the refractive 
index for the two modes, but not the mode’s dynam¬ 
ics. Simulations in the mid panels, where the thickness 
variation has been removed and the plate has a constant 
thickness /&&, show that although the So mode is not af¬ 
fected (as expected) the Ao mode it is, given that the 
field distribution is different now. 

Figure [l] lower panel, shows the normalized field dis¬ 
tribution along the central axis of the lens and for the 
four situations described before. It is clear that the 
field profile and focusing point is nearly the same for the 
two modes when the thickness of the plate is position- 
dependent, but when the thickness is kept constant the 
Ao mode is different, which clearly verifies the role of the 
thickness variation. 

The other device proposed in the present work is the 
Luneburg lens. Figure [ 2 ] shows the performance of a 
Luneburg lens of radius Rl = 10a, designed in the same 



SO mode h = h(r) 
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FIG. 1. Upper Panels: Field distribution for ka — 1 for 
a Ao polarized (left) and a So (right) polarized plane waves 
interacting with the gradient index lens where the thickness of 
the plate is designed for the simultaneous focusing of the two 
Lamb modes. Mid panels: Same as the upper panel but now 
the thickness of the plate is constant throughout the lens, so 
that the focusing of the Ao mode is different than for the So 
mode. Lower panel Normalized field distribution along the x 
axis and for y — 0 (central axis of the lens) for the So mode 
with thickness variation (black continuous line) and without 
thickness variation (blue dashed line), and for the Ao mode 
with thickness variation (red dash-doted line) and without 
thickness variation (magenta doted line). It is seen that the 
thickness variation does not affect the So mode, while for the 
Ao mode it changes the focusing point, making it identical to 
that of the So mode 


way as the GRIN lens to work for the two modes (in¬ 
clusions are now arranged in a triangular lattice, which 
is more adequate for circular clusters). The upper-left 
panel shows the real part of the normalized 2 component 
for the Ao polarization and the right panel shows the real 
part of the normalized x component for the So polariza¬ 
tion. The lower panel shows the absolute value of these 
fields along the axis of the lens. As before the wavelength 
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FIG. 2. Upper panel: Real part of the field distribution for 
the Ao (left panel) and So (right panel) polarizations inter¬ 
acting with a Luneburg lens designed by means of phononoic 
crystal with a thickness variation. Lower panel: Absolute 
value of the field distribution along the y — 0 line. Vertical 
dashed lines show the radius of the lens. 

for the two fields is ka = 1. In all the plots it can be seen 
how a focusing point near the border of the lens occurs, 
as expected. 

In the previous simulations it has been shown the 
performance of the two devices when working at the 
same wavenumber, but different frequencies, for the two 
modes. However, for a practical application it is nec¬ 
essary to see the performance of the lens at the same 
frequency. Figure [3] shows a COMSOL simulation of the 
GRIN lens working at the frequency uja/c t = 1.22, left 
panel for the Ao mode and right panel for the So one. 
We can see how the focusing point is nearly the same for 
the two modes, while it is better defined for the Ao mode 
where the wavelength is A ~ 3.6a, while for the So mode 
this wavelength is A « 7.7 a. Similar results were found 
also for the Luneburg lens. 

It must be pointed out that this type of devices is typ¬ 
ically limited by two wavelengths. In the lower limit, it 
is found that for wavelengths larger than the typical size 
of the device the diffraction dominates over refraction, 
and the device’s focusing properties are altered. In the 
upper limit, it is found that for wavelengths shorter than 
3 or 4 times the lattice constant, the phononic crystal is 
no longer homogeneous, and the effective medium model 
cannot be applied. Given that this wavenumber region 
is well defined, but that the corresponding frequency re¬ 


FIG. 3. Field distributions for the interaction of Ao mode 
(upper-left panel) and the So mode (upper-right panel) work¬ 
ing at the same frequency but different wavenumbers. Lower 
panel: Dispersion relation uj — k for the Ao (blue-continuous 
line) and the So (red-dashed line). Vertical lines show the 
wavenumber region in which the GRIN lens would work prop¬ 
erly, horizontal red-dashed and blue-continuous lines show the 
frequency region in which the device will work properly for 
the two modes. The horizontal dot-dashed line shows the 
frequency at which the simulations in the upper panels have 
been performed. 


gion is not the same for the two modes, the frequency 
region in which the device works properly for the two 
modes will be shorter than that in which it works individ¬ 
ually for each of them. The lower panel of figure [3] shows 
the dispersion relation for the So mode (red-dashed line) 
and for the A 0 mode (blue-continuous line). The two 
vertical lines define the wavenumber region in which the 
performance of the device is expected. The frequency re¬ 
gion in which the device will work properly for the two 
modes is limited below for the mode which has the larger 
wavelength, that is, the So mode, and above for the A 0 
mode, given that it has the shorter wavelength. These 
limits are shown in the figure by the red-dashed and blue- 
continuous horizontal lines, respectively. The horizontal 
black-dot-dashed line shows the frequency at which the 
simulations were performed. It is seen that there is still 
a broadband frequency region in which the device can 
work properly for the two modes. 

In summary, this work shows that simultaneous control 
of the So and Ao Lamb modes by means of a graded in¬ 
dex phononic crystal, combined with a properly designed 
thickness variation of the plate, is possible. This control 
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has been done here by designing a GRIN flat lens and 
a Luneburg lens, and it has been shown a good perfor¬ 
mance of these devices when interacting with both a So 
polarized wave and a Aq one. These devices work prop¬ 
erly in a broadband frequency region for the two modes, 
so that it is in this frequency region where most of the 
applications are expected. Other gradient index devices 
are easily envisioned, and also its extension to the control 
of the SHq mode or other type of guided modes is likely. 
This work opens the door to new type of broadband mul¬ 
timodal devices with a wide variety of applications in all 
scales. 

This work was supported by the “Agence Nationale 
de la Recherche (ANR)” and the “Delegation Generale a 
lArmement (DGA)” under the project Metactif, Grant 
No. ANR-ll-ASTR-015. Yabin Jin acknowledges 
a scholarship provided by China Scholarship Council 
No.201406260170. 


* daniel.torrent@iemn.univ-lillel .fr 

[1] M. Sarbort and T. Tyc, Journal of Optics 14, 075705 

( 2012 ). 

[2] A. A. Krokhin, J. Arriaga, and L. N. Gumen, Phys. Rev. 
Lett. 91 , 264302 (2003) 

[3] S.-C. S. Lin, T. J. Huang, J.-H. Sun, and T.-T. Wu, 
Phys. Rev. B 79 , 094302 (2009) 


[4] A. Climente, D. Torrent, and J. Sanchez-Dehesa, Ap¬ 
plied Physics Letters 97 (2010). 

[5] T. Martin, M. Nicholas, G. Orris, L.-W. Cai, D. Tor¬ 
rent, and J. Sanchez-Dehesa, Applied Physics Letters 
97 (2010). 

[6] L. Zigoneanu, B.-I. Popa, and S. A. Cummer, Physical 
Review B 84, 024305 (2011). 

[7] T. P. Martin, C. N. Layman, K. M. Moore, and G. J. 
Orris, Physical Review B 85, 161103 (2012). 

[8] S.-C. S. Lin, B. R. Tittmann, and T. J. Huang, Journal 
of applied physics 111 , 123510 (2012). 

[9] A. Climente, D. Torrent, and J. Sanchez-Dehesa, Ap¬ 
plied Physics Letters 105 , 064101 (2014). 

[10] D. Torrent, Y. Pennec, and B. Djafari-Rouhani, Journal 
of Applied Physics 116 , 224902 (2014). 

[11] A. Climente, D. Torrent, and J. Sanchez-Dehesa, Ap¬ 
plied Physics Letters 100 (2012). 

[12] O. Umnova and B. Zajamsek, in Acoustics 2012 (2012). 
[131 S. A. Cummer and D. Schurig, New Journal of Physics 

9, 45 (2007). 

[14] A. Norris and A. Shuvalov, Wave Motion 48, 525 (2011). 

[15] D. Royer and E. Dieulesaint, Elastic Waves in Solids I: 
Free and Guided Propagation , Vol. 1 (Springer Science &; 
Business Media, 2000). 

[16] T.-T. Wu, Y.-T. Chen, J.-H. Sun, S.-C. S. Lin, and T. J. 
Huang, Applied Physics Letters 98, 171911 (2011). 

[17] J. Zhao, R. Marchal, B. Bonello, and O. Boyko, Applied 
Physics Letters 101 , 261905 (2012). 

[18] X. Yan, R. Zhu, G. Huang, and F.-G. Yuan, Applied 
Physics Letters 103 , 121901 (2013). 

[19] D. Torrent, Y. Pennec, and B. Djafari-Rouhani, Physical 
Review B 90 , 104110 (2014). 



